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A STUDY  OF  THE  WAVEGUIDE  MODE 


EXCITATION  AND  PROPAGATION  IN  A 
PARALLEL  PLATE  TRANSMISSION  LINE 

ABSTRACT 

A study  of  the  waveguide  mode  excitation  and  propagation  is  made 
for  a parallel  plate  waveguide.  Since  this  configuration  is  used  in  Che 
simulation  of  the  nuclear  electromagnetic  pulse,  results  appropriate 
for  various  simulators  are  obtained  and  discussed.  The  analysis  assumes 
infinitely  wide  places,  but  Che  question  of  finite  width  plates  is 
addressed. 
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1.  Introduction 


The  parallel  plate  transmission  line  provides  a useful  means  for 
measuring  the  electromagnetic  field  Interaction.  In  this  configuration 
nearly  uniform  electric  and  magnetic  fields  are  produced,  and  Instru- 
mentation can  be  placed  beneath  one  of  the  plates.  There  Is,  however, 
a practical  problem  associated  with  the  design  and  excitation  of  the 
parallel  plate  transmission  line;  namely,  one  must  provide  a transition 
section,  typically  one  or  two  conical  transmission  line  sections,  between 
the  source  and  the  parallel  plate  configuration  (see  Figures  1 and  2). 

Insofar  as  the  transmission  line  mode  is  concerned  the  transition 
section  poses  no  problems  provided  the  conical  line  and  the  parallel 
plates  have  approximately  the  same  cross  sectional  shape  at  their  Inter- 
face and  provided  there  Is  a match  between  the  characteristic  Impedances. 
Ro%rever,  the  spherical  wavefront  propagating  along  the  conical  line  does 
not  match  the  planar  wavefront  of  the  TEM  mode  In  the  parallel  plate 
region,  hence  the  excitation  of  TE  and  TM  waveguide  modes  as  well  as  the 
TEM  mode  occurs.  Finally,  In  some  situations  a single  Input  section  Is 
not  sufficient  to  produce  the  required  field  strengths.  In  these  cases 
(usually  wide  plate  separations),  dual  Input  sections  are  used.  When 
using  dual  Input  sections  one  must  allow  provision  in  the  analysis  for 
the  fact  that  the  sources  may  not  fire  simultaneously. 

In  studying  the  interaction  of  complex  electrical  systems  with  the 
electromagnetic  pulse  generated  by  a nuclear  detonation,  the  parallel 
plate  waveguide  can  be  used  to  simulate  the  Interaction.  The  electro- 
magnetic pulse  may  be  simulated  by  driving  the  parallel  plates  from 


Schematic  Diagram  of  Two  Parallel  Plate  Simulators.  Mote  that  all 
dimensions  are  In  meters  and  the  scale  Is  1cm  > 2Sm. 


conical  Input  sections  with  appropriate  surge  generators.  There  are  a 
number  of  these  simulation  facilities  currently  in  use.  Example 
geometries  for  single  and  dual  input  parallel  plate  transmission  line 
simulators  are  shown  in  Figures  1 and  2. 

An  additional  complicating  feature  of  these  simulators  is  that  they 
are  constructed  as  finite  width  parallel  plate  transmission  lines.  The 
properties  of  the  TEM  mode  on  such  a structure  have  been  studied  exten- 
sively [1,2].  Other  studies  [3]  have  shown  that  the  waveguide  modes 

j 

for  open  region  waveguides  differ  considerably  from  the  modes  of  the 

more  familiar  closed  region  waveguides.  The  excitation  coefficients  for 

the  modes  in  an  open  region  structure  are  obtained  by  matching  boundary 

conditions  at  the  interface  between  the  conical  input  section  and  the 

parallel  plate  region  as  well  as  on  the  waveguide  walls  [3].  Apparently, 
f 

Integral  equations  must  be  solved  numerically  to  obtain  the  coefficients 
. and  propagation  constants.  Due  to  the  complexity  of  these  problems, 

the  finite-width  parallel  plate  configuration  is  approximated  by  using 
j infinitely  wide  plates.  However  expected  differences  resulting  in  the 

..  data  that  are  obtained  are  noted  and  discussed. 

* The  analysis  that  is  presented  is  validated  by  comparisons  with 

I measured  data  and  with  data  obtained  by  other  investigators. 
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2.  Electromagnetic  Fields  In  an  Infinite,  Parallel  Plate 
Waveguide  with  a Single  Conical  Input  Section 

a.  Physical  Model 

One  wishes  to  obtain  the  electric  and  magnetic  fields  In  the  parallel 
plate  region  of  the  structure  shown  In  Figure  3 In  terms  of  known  quantities 
In  the  Input  section  and  the  physical  dimensions  of  the  parallel  plates. 

Several  simplifying  assumptions  are  made.  The  plates  are  assumed  to 
be  perfectly  conducting,  the  analysis  Is  performed  In  the  plane  ASB  of 
Figure  3,  the  wavefronts  are  assumed  spherical  In  the  conical  section  and 
written  as 

E(x,z)  “ ^ 0 (1) 

where  r*  + (z  + (see  Figure  4).  Also,  It  Is  assumed  there 

are  no  reflections  In  the  conical  section  from  the  parallel  plate  region, 
and  the  electric  field  In  the  aperture  AMB,  Figure  4,  Is  assumed  to  be  a 
planar  section  of  a spherical  wavefront. 

b.  Mathematical  Model 

If  one  assumes  the  apex  of  the  conical  Input  pulser  generates  spheri- 
cal waves,  a magnetic  current  source  M - Independent  of  y-  Is  located  at 
the  plane  z ~ 0 and  M Is  given  by  [4] 

M - E X n (2) 

a 

where  E^  Is  the  electric  field  strength  at  the  aperture  or  In  the  plane  . 
z " 0.  From  the  method  of  Images, 

M - 2 E^  X n - -2(E^  • x ) i ~ y (3) 
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Figure  4.  Two-Parallel-Plate  Waveguide  with  a Single  Pulser. 
(Longitudinal  section  through  apex  S.) 


Substitution  of  (1)  Into  (3)  yields 


-jkr 


o ^ 


M - -2V  ^ 

y ° 


0 • X 


(4) 


A A 


Since  0 • X “ cos  0 = z /r  , It  follows  that 

o o 


My  " -2^0  Vo  ^ r 


(5) 


where  = x^  + z^  . The  constant  V Is  defined  by  normalizing  the 
o o o 

voltage  between  the  plates  to  unity:  i.e. 

,A 


BA 


-I 

B 


E(x,z)  • ds  = -1  V 


(6) 


Choosing  the  spherical  wavefront  through  A and  B as  the  path  of  Integration, 


l,e.,  r = £ , one  finds 
o ’ 


^jk£ 

20 


(7) 


and 


-jkC^^  + z^  - )l) 


M 


0 (x^  + z^) 
o o 


O ' V 

m 


(8) 


For  a charge  and  current  free  medium.  Maxwell's  equations  are 


7 • D - 0 


7 • B - 0 


7xE  ■ -jupH-M 


(9) 


7xH  - ju»eE 
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where  e and  u are  the  free  space  values  of  permittivity  and  permeability, 
and  the  time  dependence  is  understood  but  suppressed.  Since  M is 

y independent  and  the  plates  are  assumed  to  be  infinite  in  extent  in  the 
y direction,  Maxwell's  equations  must  be  y independent.  Thus,  for  TM  waves 


(V*  + k^)  H - jwe  M 
xz  y y 

(10) 

E - 

*^x  j 0)  e dz 

(11) 

E . _J_ 

2 j 0)  e dx 

(12) 

2 d d 

where  V = — - + — - and  k = 2Tf/A  = u)/c 
xz  dx^  dy^ 


c.  Frequency  Domain 

The  standard  technique  for  solving  (10)  is  the  Green's  function 
method;  i.e.,  is  replaced  by  G(x,2|x',0)  and  is  replaced  by 

6(x-x')  6(z).  Thus,  to  find  G , one  must  solve 


(7^  + k.*)  G(x,zlx',0)  • j la  e 6(x-x')  6(z)  (13) 

xz 

Additionally,  G must  satisfy  the  same  boundary  conditions  as  Hy  , or 


— 00 


G Is  assumed  to  be  a Fourier  series  of  even  and  odd  terms 
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I 

I 


n-0.2.4, 


I (t) 


. nTTx 

— b (z)  cos  

a I n a 


I 

n-1,3,5, , 


il*"  b°(z)  sin 
a I n a 


(14) 


with  superscripts  e and  o representing  even  and  odd  integers  respectively, 
and 

if  n ■ 0 


e - 1 

n 


if  n 0 


The  solution  of  (13)  for  the  even  terms  in  G is 


I , u)£  r ^n  mrx' 

G (x,2|x’,0)  - - — I cos 

n-0,2,4,... 


niTx  e 

cos  — 

a a 


-Jk  z 
n 


Rn 


(15) 


where  k - k^  - (nn/a)^.  To  obtain  H , one  evaluates  the  Integral 
n y 


.a/2 

Hy  (x,z)  » I My(x')  G®(x,z|x’,0)  dx' 
-a/2 


(16) 


It  follows  that 


H*  (x,z) 


z iiice 
o 


Jk£ 


®o®  n-0,2,4,... 


fr  1 ■•‘"nZ 

V fn  A®  r-s  ^ - 


(17) 


where 


_ -jk'^x'*  + z^ 

f t o 

.e  nirx'  e 

K - cos  — :i dx 

n i 


a/2 
-a/2 


x’*  + z " 


(18) 
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Similarly,  one  can  solve  for  the  odd  terms  with  the  result 


Hy  (x,z) 


z U)ee 
o 


nTTx  e 


0 a ^ 

° n-1,3,5,. 


r ,0  . nTTx 
) A sin  

“ n a 


a 

■ I 


nTTx'  e 


-jk^x'*  + z^ 


I 2 2 

X + z 


A for  single  cone  excitation  is  identically  zero;  thus 
n 


Hy(x,z) 


z b)ee 
o 


j V , . 2mrx  e 

s ) e A cos  — 

6 a 1 n n a 

o n»0,l,2,3, . . . 


In  (21)  it  should  be  observed  that  the  summation  is  now  over  all  integers, 


-jk*'x'*  + z^ 

2n'nx ' e ° , 


k*  ■ k*  ~ (2mrx/a)* 
n 


Use  of  (11)  and  (12)  in  conjunction  with  (21)  leads  to 


Ej^(x,z) 


^O®  V r ft  2mTX 

— > e A cos  e 

6 a n n a 

o n*0 


E^<x,z) 


2TTJz^e-^ 


-Jk  z 

r . . 2n'nx  e " 

) n A sin  — ; 

n»0  n 


T 
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The  first  terms  In  (21),  (24).  and  (25)  represent  the  TEM  mode  and  all 
remaining  terms  are  TM  modes.  For  finite-width  parallel  plates  TE  modes 
also  are  excited  and  propagated  [3].  Moreover  the  TE  modes  are  conside- 
rably less  attenuated  than  the  TM  modes.  For  both  sets,  the  modes  are 
attenuated  above  cut-off  due  to  the  radiation  from  the  open  sides  of 
the  waveguide.  Thus  the  results  that  are  obtained  for  the  Infinitely 
wide  plates  must  be  Interpreted  carefully. 


d.  Time  Domain.  Magnetic  Field 

Having  obtained  the  frequency  domain  results  of  section  c,  one  can 
proceed  In  several  ways  to  evaluate  the  time  domain  results.  From  (21) 
[and  also  (25)  though  this  will  not  be  discussed]  one  can  solve  the  time 


domain  problem  analytically.  From  (24),  one  can  evaluate  the  field  at 


various  frequencies  and  by  means  of  numerical  Inverse  Fourier  transform 
[5]  determine  the  behavior  in  the  time  domain. 

Showing  the  cj-dependence  expllcityly,  one  can- rewrite  (21)  as 


(26) 


(27) 


The  expression  in  (26)  is  the  frequency  response  at  a particular  a>  or 
the  response  to  a 6(t)  input.  The  Fourier  transform  of  the  step  function, 
u(t) 


is 


u(t)  - 0 


- 1 


t < 0 


t > 0 


P[u<t)l  - 


(29) 


from  which  it  follows  that  the  response  to  the  step  function  excitation 
is 


hy(x,2,t) 


w 

5 I 


u)  e du) 


(30) 


+ I h; 

n“l  ^ 


^ere  o and  n refer  to  the  zeroth  or  nth  term.  Rewriting  the  h^  and 
collecting  terms  in  a convenient  manner,  one  obtains 


h“° 

y 


dx’ 

(31) 


and 

f(t,t,z',x')  - t + (i  - z)/c  - *x'*  + z*/c 

Thus,  (31)  becomes 


1 

o! 

M 

fa/2 

1 

/x'*  + z* 

J 

u T •• 

c c 

, 4 

o 


dx' 


(32) 
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and,  since  I ^ /x'*  + z*  , then  if  t ^ z/c 


so 


1 

a n 


t > z/c 


(33) 


In  general  one  must  evaluate  (32)  subject  to  the  condition  that 


ct  > '^x'*  + z * + z - I 

— o 


Since  x*  varies,  this  means  the  Integral  must  be  done  in  one  of  three 
ways  depending  on  the  time,  t . 

1.  If  ct  < z^  + z - 2.  , the  Integral  in  (32)  is  identically 


zero. 


2.  If  z + z - £ < ct  < -/I  « 
o W 1 2 


4 


+ z^^  + z - £ , the  upper  limit  on 


this  integral  is 

X - ‘'(ct  + £ - z)'  - 


If  ct  — ^ 
is  a/2  . 


f? — i 

' a 

+ z^ 

l2j 

o 

These  statements  are  basically  casuallty  conditions. 

Using  the  Identity  6(a))  (Ko))  » 6(a))  (()(0)  , one  finds 


sn 


ra/c  2n'n’x' 

I cos  


dx'  • — 

,£.2  + ^2  2Trj 

o 


r 


x'*  + z* 
o 


VfWI 


(34) 


d(^)  . 


» 


■ I 


The  second  integral  In  (34)  can  be  transformed  analytically  [6]  using 
the  Identify 


1 

2irj 


- y* 
- y^ 


du)  • J (y^t^  - X*  ) , t > 


(35) 


The  final  result  of  these  manipulations  is 


O9  ^ 

1 ^ o V 2nTrx 

^ o n*l 


^x. 


cos 


2nxx 


• 1 


a x'"  + zl 


J(ct  + i - /^^  + zj  )^  - ^ dx- 


(36) 


Again,  the  Integral  in  (36)  must  be  evaluated  subject  to  the  same  con- 
ditions discussed  for  (32),  and  f^  depends  on  the  evaluation  of  (32) 
(e.g.,  if  ct  ^ /[f  + z - ^ . then  fi  - 1 .).  The  expression 

for  hy  , (36),  must  be  evaluated  numerically.  First,  define  a new  func- 

tlon  'P(N)  which  is  identical  to  h^  except  the  summation  runs  from 
one  to  N rather  than  from  one  to  infinity.  The  function  't'(N)  versus  N 
is  shown  in  Figure  5.  The  function  calculated  in  the  program  is  defined 


as 


40 

I 'i'(N) 

.8  , ..  N-20 

h (x,z,t)  ■ 

y 21 


(37) 


One  other  comnent  is  appropriate.  The  Integrals  in  (37)  are  evaluated 
using  a 15-polnt  Gauss-Legendre  quadrature  formula.  As  n increases, 
the  function  to  be  evaluated  oscillates  more  rapidly.  Therefore,  a change 
of  variable  is  performed  such  that  as  n Increases,  the  number  of  Gauss- 
Lengendre  zones  from  x'*0  to  x'*a/2  Increases  proportional  to  n. 
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The  y component  of  the  magnetic  field,  h®(x.z.t),  doe  to  a step 
function  excitation  Is  plotted  with  respect  to  time  (t)  and  shown  In 
Figures  6 for  several  field  points  In  an  ALECS-type  facility,  while 
Figure  7 Is  appropriate  for  the  ARES  geometry.  Note  that  time  Is 
measured  from  t - 0 . the  time  at  which  all  the  dipoles  In  the  plane 
2-0  have  turned  on.  Corresponding  frequency  domain  (Fourier  Trans- 
forms of  hy)  are  shown  In  Figures  8 and  9.  Note  that  singularities  occur 
at  the  cavity  resonances  of  the  parallel  plates. 


Time  Domain.  Electric  Field 

The  electric  field  calculations  in  the  time  domain  have  been  per- 
formed using  the  Inverse  Fourier  transform  rather  than  an  analytic  solu- 
tion. E^(x,2,U)),  (24),  has  been  evaluated  for  100  different  frequencies 

varying  from  the  d-e  term  to  100  MHz.  The  frequencies  chosen  were  as 
follows: 


f 

Af 

0 

- 

.5  - 1 MHz 

.1  MHz 

1 - 10  MHz 

.2  MHz 

10  - 100  MHz 

2 MHz 

e,(x.z.t)  has  been  evaluated  using  standard  Fourier  transform  techniques: 
In  this  procedure  It  Is  assumed  that  E^fx.z.co)  Is  a linear  function  of 
frequency  between  the  Individual  calculated  field  values  [5].  Since 
e^(x,z,u>)  as  given  by  (24)  Is  the  Impulse  response  then  It  must  be 
multiplied  by  the  Fourier  transform  of  the  actual  pulse  driving  the 
plates-ln  the  previous  section  this  pulse  was  assumed  to  be  a unit  step. 
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Figure  6:  The  horizontal  magnetic  field  produced  by  a unit  step 
excitation  of  the  ALECS  simulator. 


. 10 


Figure  8;  The  Fourier  transform  of  the  horizontal  magnetic 
field  produced  by  a unit  step  excitation  of  the 
AI.FCS  simulator. 


f in  ftHz 


Figure  9 


The  Fourier  tr/insform  of  the  horizontal  magnetic 
field  produced  by  a unit  step  excitation  of  the 
ARES  simulator. 


I 


field  produced  by  a unit  step  (with  exponential 
decay)  excitation  of  the  ALECS  simulator. 
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Figure  12;  The  ' ourler  transiorm  of  the  vertical  electric 
flelc  produced  by  a unit  step  (with  exponential 
decay)  excitation  of  the  ALECS  simulator. 
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Figure  13:  The  Fourier  transform  of  the  vertical  electric 
J field  produced  by  a unit  step  (with  exponential 

I decay)  excitation  of  the  ARES  simulator. 
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Figure  14:  The  Fourier  transform  of  the  vertical  electric 
field  produced  by  a unit  step  (with  exponential 
decay)  excitation  of  the  ARES  simulator. 
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The  vertical  electric  field  produced  by  a unit  step  (with 
exponential  decay)  excitation  of  the  ALECS  simulator. 


Figure  17;  The  vertical  electric  field  produced  by  a unit  step  (with 
exponential  decay)  excitation  of  the  ALECS  simulator. 


Figure  18:  The  vertical  electric  field  produced  by  a unit  step  (with 
exponential  decay)  excitation  of  the  ALECS  simulator. 


The  vertical  electric  field  produced  by  a unit  step  (with 
exponential  decay)  excitation  of  the  ARES  simulator. 


f ( ii.  MHz) 


Figure  23:  Evaluation  of  the  CW  performance  of  the  ALECS 
facility-magnetic  field  determination  at  the 
center  of  the  working  volume,  x=0,  z=7.5. 


However,  for  convenience  in  applying  the  inverse  Fourier  transform  here 
the  unit  step  with  exponential  damping  is  considered,  i.e.,  u(t)e 
where  a*  4.8  x 10*  s 

The  complex  magnitude  of  the  Fourier  transform  of  the  vertical 
electric  field  is  shown  in  Figures  10  through  15.  In  contrast  with  the 
results  of  Figures  8 and  9 the  Fourier  transform  of  the  vertical  electric 
field  does  not  exhibit  slngu’ arities  at  the  cavity  resonances.  Results 
appropriate  for  the  ALECS  and  ARES  simulators  are  given.  Note  that  it 
is  only  when  the  frequency  is  above  the  lowest  cut-off  frequency  that 
there  is  any  z-dependence  of  the  fields. 

Time  domain  results  for  the  vertical  electric  field  are  shown  in 
Figures  16  through  21.  These  indicate  that  the  TEM  mode  contribution 
is  dominant. 

f . Comparison  of  Experimental  Data  and  Theoretical  Results 

for  the  ALECS  Facility 

In  1974  Giles  et.al.  [7]  evaluated  the  performance  of  the  ALECS 
facility  in  the  CW  mode.  In  particular,  B and  D measurements  were  made 
versus  frequency  at  the  center  of  the  working  volume.  Figures  22  and  23 
shown  the  results  of  these  measurements.  On  the  same  figures  one  can 
see  the  results  of  the  theoretical  analysis  for  CW  performance.  In 
order  to  compare  results,  and  were  chosen  such  that  the 

experimental  and  theoretical  results  matched  at  5 MHz. 

Before  Interpreting  the  foregoing  results,  the  question  of  finite 
width  versus  infinite  width  plates  should  be  considered.  Marin  [8] 
found  that  for  finite  width  plates  the  TE  modes  were  most  important, 
i.e.,  they  exhibited  less  attenuation  above  cutoff,  and  the  analysis 


used  here  did  not  include  TE  modes.  However,  to  propagate  TE  modes  the 
plates  must  support  transverse  currents  and  since  the  ALECS  simulator 
(as  well  as  the  ARES  simulator)  is  constructed  with  parallel  wires 
rather  than  actual  plates,  the  ALECS  simulator  will  not  propagate  TE 
modes.  And  also  since  the  TEM  mode  contributes  very  significantly  to 
the  simulator  field  there  is  good  agreement  between  the  measured  and 
calculated  data.  The  disparity  increases  significantly  at  higher 
frequencies  primarily  because  the  waveguide  modes  supported  by  the 
actual  finite  width  waveguide  are  evanescent.  The  infinities  in  the 
calculated  results  occur  since  there  are  assumed  to  be  no  reflections 
in  the  conical  input  section. 
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3.  Electromagnetic  Fields  in  an  Infinite,  Parallel 


Plate  Waveguide  with  a Double  Conical  Input  Section 


a.  Physical  Model 

In  some  Instances  one  wishes  to  obtain  field  values  in  the  parallel 
plate  region  of  a simulator  that  are  higher  than  those  attainable  from  a 
single  pulser.  In  these  situations,  one  solution  is  to  use  two  separate 
pulsers  and  two  input  sections.  Figure  24.  One  problem  that  arises  is 
the  nonsimultaneous  firing  of  the  sources,  and  the  analysis  must  take  this 
into  account  if  it  is  to  provide  meaningful  data. 


b.  Mathematical  Model 

As  in  section  2,  each  source  is  assumed  to  generate  a spherical 
wave , E 


E = V 


o r 


where  r = (x  + a/4)  + (z  + z^)  for  the  upper  (-)  and  lower  (+)  cones. 

Again,  normalizing  each  source  to  1 volt,  one  obtains  a magnetic  source 


for  the  upper  cone  as 


-jk  - a/4)^  + zj  - 

0 [ (x  - a /4 ) ^ + z ^ ] 
o o 


For  the  lower  cone,  the  magnetic  source  becomes 


-Jk  |'^(x  + a/4)*  + z^*  - X-j 


0 l(x  + a/4)^  + z M 
o o 
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Figure  24;  A Design  for  a Parallel-Plate  Transmission  Line 

Simulators.  (View  parallel  to  plates  and  perpendicular 
to  direction  of  propagation,  or  z-direction. ) 
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Substitution  of  (39)  into  (16)  [the  upper  source  analysis  is  performed 


independently  of  the  lower  source]  leads  to 


“H®(x,z,a)) 


z eu 

_2_  y 

0 a n 'I  ! 

\ n-0,2,4,., 


E 
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, z > 0 


where 
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a 
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-jk  (x'  - a/4)^  + z‘ 
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, z > 0 
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-jk  *'  (x'  - a/ A)^  + z 


° dx’ 


(x’  -a/A)^  + z^^ 


0 o 

For  the  lower  cone  one  obtains  expressions  for  and  that 

are  Identical  to  (Al)  and  (A3)  except  and  A^  must  be  replaced 

by  A and  A where 
^ n n 


-jk  (x'+  a/A)^  + 

o dx' 


(x' + a/A)'  + z] 


-jk  (x'  + a/A)  * + z 


o dx' 


(x'  + a/A)*  + z^ 
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For  nonslmultaneous  firing  of  the  sources,  (41)  and  (43)  should  be  modi- 
fied so  that  they  differ  in  phase  by  OJT  where  T is  the  time  duration 
between  the  firings, 

c.  Time  Domain.  Magnetic  Field 

The  step  response  for  the  even  terms  in  the  upper  cone  excitation 
is  obtained  as  follows.  Expression  (41)  is  multiplied  by  the  frequency 
spectrum  of  the  step  function 


TT6(a))  + — 


and  the  Fourier  transform  is 


u^se 
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2ir 
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(47) 


As  in  the  single  cone  case  in  section  2,  the  o)  Integral  can  be  performed 
analytically  leading  to 


“hy®(x,z,t) 


ju  + _::£ 

an  a6 
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o n“2,4,. 


cos 


nTTx 


i/2 


cos 
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+ ^ - »'(x’  -a/4)*  + z 
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(48) 
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with  the  restriction  that  ct  + i (x*  -a/4)^  + > z . In  a 

o — 

similar  fashion  odd  terms  for  the  upper  source,  ^h®^  , and  even  and  odd 

IL  se  A so 

terms  for  the  lower  source,  and  respectively,  can  be  derived. 

Finally,  if  one  allows  the  upper  source  to  fire  after  the  lower  source  by 
a time  interval  T , one  obtains 


, s,  u.  se.  . u,  so(x,z,t) 

hy(x,z,t)  - hy  (x,z,t)  + hy  ' ’ ’ ' 

, se.  , . 2.,  so,  X 

+ h (x,z,Ti)+  h (x,z,Tj) 


where  Tj  = t + T . In  detail.  (49)  becomes 
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The  coiranents  regarding  the  solution  of  (36)  in  section  2 are  also 
applicable  to  the  solution  of  (50). 


d.  Time  Domain,  Electric  Field 

u s s 

Having  derived  H^(x,z,<i))  and  Hy(x,z,a)),  one  can  obtain  the 
Fourier  transform  of  the  vertical  electric  field  E^(x,z,a))  from  (24). 
Then,  in  exactly  the  same  procedure  as  that  described  in  section  2e, 
time  domain  electric  fields  may  be  obtained  using  the  Inverse  Fourier 
transform.  Also  as  before  the  electric  field  calculation  is  made  for  an 
exponentially-decaying,  unit-step  excitation. 


e.  Numerical  Results  for  the  ATLAS  I Simulator 

While  the  ALECS  and  the  ARES  simulators  could  not  propagate  TE  modes, 
the  ATLAS  I simulator  can  since  the  waveguide  plates  are  a wire  mesh  and 
hence  can  support  currents  transverse  to  the  direction  of  propagation. 
Therefore,  the  results  to  be  presented  subsequently  are  valid  only  in  a 
qualitative  sense  since  the  analysis  did  not  include  TE  modes.  However, 
the  excitation  of  the  TEM  mode  and  the  asynchronus  firing  of  the  pulses 
should  be  treated  quite  accurately  by  the  analysis.  Accordingly  a 
copious  amount  of  data  is  presented.  Since  the  ATI.AS  I simulator  is  still 
under  construction  there  is  no  measured  data  to  validate  the  analysis. 

For  a detailed  discussion  of  the  predicted  behavior  of  the  ATLAS  I 
simulator  the  authors  defer  to  references  [9]  through  [12]. 

First  the  Fourier  transform  of  the  vertical  electric  field  is  shown 
for  various  positions  within  the  ATLAS  I simulator  with  an  assumed 


typical  asynchronus  firing  of  the  pulsers.  Magnitude  plots  are  given  in 


Figures  25  through  28  for  various  points  within  the  working  volume. 


r 
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Figure  25:  The  Fourier  transform  of  the  vertical  electric  field 
produced  by  a unit  step  (with  exponential  decay) 
excitation  of  the  ATLAS  I simulator. 


Figure  26:  The  Fourier  transform  of  the  vertical  electric 
field  produced  by  a unit  step  (with  exponential 
decay)  excitation  of  the  ATLAS  I simulator. 


(x,z,a))I  In  nV/m*Hz 


Figure  29:  The  Fourier  transfcin  of  the  vertical  electric 
field  produced  by  a unit  step  (with  exponential 
decay)  excitation  of  the  ATLAS  I simulator. 


In  mV/ci 


Figure  31:  Tlie  vertical  electric  field  produced  by  a unit  step 
(with  exponential  decay)  excitation  of  the  ATLAS  1 
simulator. 


For  a unit  step  excitation  (even  with  exponential  decay),  the  waveguide 
inodes  are  much  more  significant  for  the  ATLAS  simulator  than  for  the 
ALECS  or  ARES  due  to  its  dual  input  sections  and  due  to  its  lowest  cut- 
off frequency  for  the  waveguide  modes  being  only  5.7  MHz.  It  is  readily 
noted  that  the  higher  order  modes  become  increasingly  significant  with 
distance  from  the  input  transition  section.  If  only  a pure  TEM  mode  were 
excited  in  the  simulator  then  the  magnitude  plots  of  the  Fourier  trans- 
form would  be  straight  lines. 

A significant  contribution  to  the  excitation  of  the  higher  order 
modes  is  the  asynchronous  firing  of  the  pulsers  driving  the  input  sections 
of  the  ATLAS  simulator.  This  can  be  seen  from  figure  29.  However,  the 
deep  nulls  in  the  curves  at  50  MHz  and  100  MHz  are  due  to  dlstructive 
Interference  occurring  as  a result  of  the  phase  difference  between  the 
two  sources  of  the  two  input  sections 

As  mentioned  in  the  foregoing  the  ATLAS  I simulator  can  support  TE 
modes  whereas  the  analysis  includes  only  TM  modes.  For  finite  width 
plates  both  the  TE  and  TM  modes  are  evanescent,  but  for  infinite  plates 
the  modes  do  not  radiate.  Hence  the  effect  of  the  TM  modes  is  over 
emphasized  by  the  analysis  (as  shown  in  Figures  22  and  23)  and  the  effect 
of  the  TE  modes  under  emphasized.  One  would  expect  then  that  the  rela- 
tive Importance  of  the  higher  order  modes  as  exhibited  in  Figures  25 
through  28  should  be  representative.  But  the  asynchronous  effect  as 
shown  in  Figure  29  should  be  quite  accurate. 

In  Figures  30  through  33  the  time  domain  vertical  electric  field 
Is  presented  for  an  exponentially  decaying  unit  step  excitation  of  the 
dual  input  sections  to  the  ATLAS  I simulator.  Note  that  at  the  center 
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Figure  34:  The  vertical  field  produced  by  a unit  step  (with 
exponential  decay)  excitation  of  the  ATLAS  I 
simulator. 


/ 


of  the  working  volume  (x=0,z=50)  the  peak  electric  field  is  about 
26  sllllvolts/raeter . Baum,  et  al.  [12]  predict  that  the  value  should 
be  about  20  rallllvolts/meter . In  arriving  at  their  result  Baum,  et  al. 
use  ray  theory  primarily.  It  is  further  noted  that  the  overall  time 
domain  response  follows  the  prediction  presented  by  Baum,  et  al,  con- 
sidering that  their  assumed  excitation  is  a unit  step.  Also  the  general 
behavior  at  x = 13,  z = 50  follows  the  roughly  sketched  waveform  that 
Baum,  et  al.  predict  for  their  point  c.  They  predict  an  initial  peak 
of  about  11  mllllvolts/meter  whereas  the  present  analysis  yields  about 
17  millivolts/meter,  which  is  quite  satisfactory  agreement  in  view  of 
the  completely  different  analyses  that  are  used. 

In  their  predictions  of  the  electric  field  produced  by  the  ATLAS  I 
simulator  Baum,  et  al . [12]  do  not  consider  asynchronous  firing  of  the 
Input  pulsers.  The  effect  in  the  electric  field  at  the  center  of  the 
working  volume  due  to  the  nonsimultaneity  of  the  excitations  of  the  dual 
input  sections  is  exhibited  in  Figure  34.  First  there  Is  an  Increase 
in  the  rise  time  and  second  if  the  asynchronous  time  is  sufficiently 
large,  e.g.  T = 10  ns  , a precusor  peak  occurs.  These  precusor  peaks 
are  even  more  pronounced  in  the  horizontal  magnetic  field  computations 
as  shown  in  Figure  35  for  the  same  location  In  the  simulator.  The  magni- 
tude of  the  Fourier  transform  of  the  horizontal  magnetic  field  is 
shown  in  Figure  36  for  a typical  location  in  the  ATIJ^S  I simulator. 
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Figure  35;  The  horizontal  magnetic  field  produced  by  a unit 
step  excitation  of  the  ATLAS  I simulator. 
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